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Abstract 

Motivated by the concept of shape invariance in supersymmetric quantum mechan- 
ics, we obtain potentials whose spectrum consists of two shifted sets of equally spaced 
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the singular term oix~'^ always falls in the transition region —\<a.<\ and there is a 
(5-function singularity at a; = 0. 
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It is well-known that a set of equally spaced energy eigenvalues corresponds to a particle 
moving in a harmonic oscillator potential. More specifically, in units of ?i = 2m = 1, the 
set of energy levels 

En=nuj (n = 0,l,2,...) (1) 
corresponds to the one dimensional potential 

V{x) = ^w^x^ - ^oj (-0O < X <oo) . (2) 

However, if the particle is not allowed to move along the entire j;-axis and is constrained to 
remain on the half-axis, a simple potential with eigenvalues given by eq. (||) is 



— LO^r^ + ^ ^ ' / + - (0<r<oo; l> 

16 r2 2 V 2/ ^ ' - 2' 



These potentials are not unique. It is well-known that they can be systematically deformed 
to produce multi-parameter isospectral families either via inverse scattering methods or 
using supersymmetric quantum mechanics techniques 

Now, instead of a single set of equally spaced eigenvalues, let us make an innocuous 
extension and ask for potentials which possess two shifted sets of equally spaced eigenvalues. 
The energy levels are given by n(a;o -|- toi) and n{ujQ + oJi) +ujq with n = 0, 1, 2, • • • . Clearly, 
there are two sets of energy levels with equal spacing {ujq + liJi), with one set displaced by 
(Jo with respect to the other set [see Fig. 1]. To our knowledge, potentials which possess the 
spectrum of Fig. 1 have not been discussed in the literature. We will obtain simple, explicit 
potentials using an approach motivated by recent developments |^ ^ in obtaining new 
shape invariant potentials in supersymmetric quantum mechanics |Q, ^. These advances 
have been made using novel choices of the function / appearing in the change of parameters 
ai = f{ao) in the shape invariance condition. In this paper, we use a transformation 
such that /^(ao) = oq. The resulting shape invariant potentials are found to be singular 
but analytically solvable and will have a spectrum consisting of two shifted sets of equally 
spaced energy levels. Other potentials with the same spectrum can be readily obtained by 
making isospectral deformations [Q, |^], but this aspect will not be further developed here. 
We discuss applications to multiparticle problems |9|. In particular, we find that there 
exist new eigenstates interlaced between the usual equally spaced Calogero spectrum of 
eigenvalues. These new states are a direct consequence of the fact that potentials with a 
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OLX ^ singularity with strength a in the transition region — | < a < |, possess additional 
bound states with singular but normalizable wave functions. 



Shape Invariant Potentials: 

Let us recall that supersymmetric partner potentials V±(x,ao) are related to the super- 
potential PF(x,ao) by 

y±(x,ao) = H^'(x,ao)±PF'(x,ao) , (4) 

where oq is a set of parameters. These partner potentials are shape invariant if they both 
have the same x-dependence upto a change of parameters a\ = /(ao) and an additive 
constant R{aQ). The shape invariance condition is 

V+{x,ao) = V-{x,ai) + R(ao) , (5) 

or equivalently 

W'^{x,ao) + W'{x,ao) = W'^{x,ai)-W'{x,ai) + R{ao) . (6) 
The property of shape invariance permits an immediate analytic determination of energy 



eigenvalues eigenfunctions [ pX| ] and scattering matrices [11 1. For unbroken supersymme- 
try, the eigenstates of the potential V- (x) are: 

n-l 



ak 



k=0 



oc e , ao) = [-^ + ao)] «i) = 1, 2, 3, ...) . (7) 

Note the existence of a zero energy ground state, characteristic of unbroken super symmetry. 



Change of Parameters: 

Until 1992, the only known solutions of the shape invariance condition corresponded to 
a translational change of parameters 

ai = /(ao) = ao + (3 , a„ = /"(ao) = ao + n/3 . 

Subsequently, the scaling change of parameters 

ai = /(ao) = qao , a„ = /"(ao) = q'^ao 



was studied and shown to produce new classes of shape invariant potentials 1^, ^, which 
included self-similar potentials ||l^ as a special case. The above mentioned change of param- 
eters are shown in Fig. 2. In general, the changes are monotonic, and the parameters never 
acquire the same values again. The only exceptions are the special cases /3 = and q = 1, 
both of which have /(ao) = clq and hence a„ = On-i = • • • = oi = oq. For this situation, 
the energy levels are equally spaced since ^ — -E^Ji = -R(ao)- If one denotes the constant 
value of R{ao) by oj, then the shape invariance condition [eq. (^)] reads W'{x,ao) = u)/2. 
This gives W{x,ao) = ujx/2 and the standard harmonic oscillator potential y_(x) of eq. 

Let us now consider a new change of parameters defined by /^(ao) = clq- Consequently, 
ao = 02 = ^4 = • • • , /(ao) = ai = 03 = as = • • • . 

There are many functions which meet the above requirement. One such function is /(ao) = 
C/ao, where C is any real constant. Another choice is /(oq) = 2^ log[coth(Cao)]. The 
energy levels will be spaced alternately by i?(ao) = ujq and R{ai) = uji. This spectrum 
corresponds to two shifted sets of equally spaced eigenvalues. 

Single Particle Potentials: 

To explicitly compute the potentials associated with two shifted sets of equally spaced 
eigenvalues, one needs to solve the shape invariance conditions: 

W\x, ao) + W'{x, ao) = W^{x, ai) - W'{x, ai) + too , 

W^{x,ai) + W'{x,ai) = W'^{x,ao) - W'{x,ao) +iOi . (8) 

The solution is readily obtained by straightforward manipulations. The superpotential 
W{x,ao) is given by 

{lvq + u;i)x {u;o - uji) , . 

^("'"°) = 4 + 2x(a;o + a;0 ' 

and the corresponding potential V-^ (x) is 

/ X (wo - wi)(3tJo + ^^1) t^i , (a;o + wi)^a;2 
V-(x,ao) = ^ . (10) 

^ ^ 4(wo+^i)V 2 16 ^ ^ 

The above results are correct everywhere except at the point x = 0, where the superpo- 
tential W{x, ao) has an infinite discontinuity. Such a discontinuity is not acceptable, since 
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while writing eqs. @ to (^), we have implicitly assumed the existence of a continuous 
superpotential with a well-defined derivative giving rise to unbroken supersymmetry and a 
zero energy ground state. To study the problem further, consider a regularised, continuous 
superpotential W{x, uq, e) which reduces to W{x, uq) in the limit e ^ 0. One such choice is 

W{x,ao,e) = W{x,ao) f{x,e) (11) 

where 

f{x,e) = tanh^ ^ . (12) 

The moderating factor / provides a smooth interpolation through the discontinuity, since 
it is unity everywhere except in a small region of order e around x = 0. In this region, 
W{x,ao,€) is linear with a slope 2t'^{ujo'+'^i) '^^^'^^ becomes larger as e gets smaller. The 
potential VL(x,ao,e) corresponding to the superpotential W{x,aQ,e) is 

V-{x,ao,e)=W\x,ao,e)-W'{x,ao,e) . (13) 

In the limit e ^ 0, V-{x, oq, e) reduces to 

V. (x, ao) = V. (x, ao) - 4 W{x, ao) [20(x) - 1] 5{x) , (14) 

where 9{x) is the unit step function and we have used lim^^o ^sech^| = 5{x), lime^o tanh | = 
[29{x) — 1]. Thus we see that the potential V-{x,ao) has an an additional singularity at 
the origin over y_(x,ao) given by ^}{x) = — [26{x) — 1] 5{x). Under scaling this 

singularity transforms as Q{Xx) = ■^Q{x) (A > 0). The sign of this singularity depends 
upon the values of ujo and wi. It is best to consider the two cases loq > uji and loq < uji 
separately. 

Case 1: ujq > uji. An example of a superpotential for this case is shown in Fig. 3(a). The 
corresponding potential is drawn in Fig. 1(a), along with a few low lying energy levels. Note 
the sharp attractive shape of the potential near x = 0. In the limit e — > 0, this attractive 
(5-function singularity at the origin is instrumental in producing a bound state at £"0 = 0. 
Case 2: loq < uji. An example of a superpotential for this case is shown in Fig. 3(b). Again, 
proceeding as in case 1, one gets a repulsive singularity in the potential at the origin. This 
is shown in Fig. 1(b). 

Naively, in the limit e ^ 0, the potential of eq. ([To| ) appears identical to a three 
dimensional oscillator with a frequency to = ^'^o+^i) g^j^^ angular momentum / given by 
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l{l + 1) = ^^°4^^^i'_|j^'^°^^^'' • However, there are some subtle but important differences. It 
is defined over the entire real axis (— cxd < x < oo) and not just the half line, and a very 
interesting feature of this potential is the inevitable "softness" of the inverse square term. 
For all values of parameters loq and uji , the coefficient of the inverse square term is bounded 
between —j and |; thus the potential falls exactly in what is known as the transition 
region [l5| . More specifically, for case 1 {ojo > wi), one has < + 1) < | and 
for case 2 [ujq < "^i) one has —j < l{l + 1) < 0. The important special case of the one 
dimensional harmonic oscillator has ujq = uJi: it corresponds to /(/ + 1) = and no 
singularity. For transition potentials, both solutions of the Schrodinger equation are square 
integrable at the origin. If one seeks eigenstates in a semi-infinite domain (0 < x < oo), one 
just keeps the less "singular" solution of the two[13, 14, 15]. However, for a fully infinite 



domain (— oo < x < oo) with a soft singularity at the origin, both are acceptable square 
integrable solutions of the Schrodinger equation, and must be retained to form a complete 
set. Single particle eigenstates for the potential y_(x,ao) = y_(x,ao) + ^(x), where il(x) 
represents the point singularity at the origin, are readily obtained from eq. The lowest 
four are: 



Eq = Q; Vo oc X 2{-o+"Ji) g-Ka^o+t^l)^' ^ (15) 

El = ujq; Vi oc x^+2{?o+-i) g-K^o+'^i)^' , 

E, = uo+u:i; 02 oc ( - 1 + ^^R±^A ^-^^ , 

i?3 = 2^o + a;i; 03 oc f-^^l^ - 3 + ^^^^^x^) x^+5^ e-|(-o+-)-^ 

In Fig. 4, we have drawn the lowest few eigenfunctions for two choices of the parameters 
Wo and uji. Note that for ujq > cji, the ground state eigenfunction diverges at the origin, 
whereas for loq < wi it vanishes. For the intermediate situation loq = loi, the ground 
state is finite at the origin and is just the standard Gaussian solution of a one dimensional 
oscillator. In Fig. 5, we have plotted the five lowest energy eigenvalues Eq to E^ as a 
function of a = {loq — loi){2>ujq + loi)/^{loq + ui)^ holding lo = {loq + uji)/2 fixed, a is the 
coefficient of the singular term in eq. (|^). 2a; is the equal energy spacing in the two sets of 
energy levels. General expressions for these eigenfunctions and corresponding eigenenergies 
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are given by 



E2n = n{ljjQ + UJi) , 
E2n+1 = n{uJQ + OJi) + UJQ 



UJX 



UJX 



where d 



to 



, and Ln are the standard Laguerre polynomials. 



(16) 



Two Particle Case: 

Consider the case of two particles interacting via the potential VI (^J^^, qq) of eq. (p^, 
consisting of the Calogero part V-{^^^£^, oq) and the additional 5- function singularity. The 
Calogero part [eq. (10)] has been extensively discussed in the context of solvable many 
body problems in one dimension However, as we mentioned earlier, if the coefficient 
a of the singular term in eq. ( |lO|) is chosen to be between — ^ and |, the "softness" of 
the singularity permits the penetration of particles through the barrier at the origin. For 
the two particle case this means that particles are able to interchange their positions by 
passing through each other. This leads to the appearance of additional eigenstates that are 
absent for a potential with a hard singularity (a > |), where only those wave functions 
that vanish at the origin are allowed. As a consequence of this "tunneling", eigenstates 
break into groups of two, with the splitting of the two states vanishing as the coefficient 
approaches |. However, in our case, for the potential VI ( ^^^^ , oq), particles also see the 
point singularity 17( ^'^^^ ) as they pass through each other. The corresponding Schrodinger 
equation is 



52 



dxi dx\ 



+ 



{Xi - X2) + 



2a 



(xi - X2Y 



Xl - X2 

~72~ 



%Ij{xi,X2) = Elp{xi,X2). 



[17) 



Changing variables to center of mass coordinate R = and relative coordinate x 

and then eliminating the motion of the center of mass, we arrive at 

v2 



(18) 



The domain of x is the entire real line. Thus, as discussed earlier for the single particle 
case, both solutions of the Schrodinger equation are necessary to generate a complete set 
of eigenstates on the full line. Now with the identification loq = 2auj and uJi = 2(1 — a)u;. 



where a = + 4a, eq. (|l8|) reduces essentially to eq. (^). Prom eq. (16), the energy 
eigenvalues are given by 



E- 



2n 



2nuj, E2n+i = 'inuj + 2abJ . 
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Three Particle Case: 

The Schrodinger equation for a three particle system that we consider here is: 

) 

where the singular function [x) is given by Q [x) = _M2£Mi:il (^(x), where h is defined by 
a = b{b + 1). To avoid "falling to the center" [|l3| , we shall confine our attention to a > — ;|. 
Transforming to center of mass and relative coordinates 

X1+X2+X3 Xi- X2 X1+X2- 2X3 

R= ^ ,x=—^,y= , (20) 

and then eliminating the center of mass variable, we get 
f 92 92 3 w <, o\ / 1 4 



Variables x and y span a two dimensional configurational plane. We parametrize this plane 



in polar coordinates r = ^/x^ + ]p and = tan ^ |. The range of these new variables are 
given by < r < 00 and < (p < 2tt. The Schrodinger equation in the variables r and (j) is 

f d'^ 19 822 if 5^ 




where we have used Q (j-sm{(j) + ^^tt)^ = ^sin((/> + ^vr)^. Owing to the zeros of the 
function sin((/) + ^vr), the fi-singularities are at ^ (/c = 0, 1, • • • , 5) and they can be 
written as J2k=o^{4' ~ Substituting \l:{r,(f)) = R{r)f{<p), and then using the usual 

procedure for separation of variables, we get: 

^ ^--^ + -yr' + f]Rir)=mr), (23) 



r dr 8 r 



and 



f 1 1 1 \ 



sm sm^(0+|7rj sm^((p+|7r) 

n=2 / n ^ 

2n 



n=0 



+ ^l^sin(</) + -7r) /,(</>) = 6F/^(0). 



The radial equation is same as that of a three dimensional oscillator. Its eigenfunctions and 
eigenvalues are given by 



Rnl{r) = r'^exp 



1 /3 2 



Lb 



—u)r 



En = ^-{2n + bi + l),Lo (25) 



where are Laguerre polynomials. To determine energy eigenvalues En, we still have to 
solve eq. (pi). We first make a simplification using the identity [0] 



+ 



+ 
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sm 



sm 



sin^((/) + Ivr) 



sin^(3(/)) 



Thus eq. (pi) now reduces to 



9a 



n=5 



k2 + „;„2 



■ ,„,+9E^(3'A-nvr) /,(,/.) =62 /^(<^) 



(26) 



where once again we have used the scaling property of the singular function fi. Since the 
effective potential of eq. ([2^) is periodic with a period vr/S, we only need to obtain a 
solution in the domain < (/> < ^. Using even and oddness of eigenfunctions, solutions can 
be extended beyond </) = 7r/3, as lucidly discussed in ref. [p. Changing variables to ^ = 3(/) 
yields 

+ + (0 + (e - ^)| fm = Bf MO, where = . (27) 

The new domain is < ^ < vr. The major role played by the fi-singularity at the origin 
is to properly describe the discontinuity suffered by derivatives of the all even-states at 
^ = and ^ = vr. Now we perform another change of the independent variable from ^ to 
z = cos^ |, and define fi{z) = [z{l — z)]^ h{z), where 5 will be chosen judiciously to remove 



the singularity in eq. (|27|). For 5 



i±4^,weget 



25 + - - (1 + 45)z 



dh 
dz 



46^ - B? 



h{z) = 0. 



(28) 



This has the form of a hypergeometric equation and the general solution is 



+ C2Z 



ciz\1-zYf 
i-5 



25 + Bu25 - Bu- + 25,z 



[l-zfF 



\ + Bu\-Bul-26,z 



(29) 



Usually one of these solutions has an unacceptable singularity near z ~ 0, and the other 
solution is retained. However, for softly singular cases {-\ < a < |) both solutions are 
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normalizable, and necessary to generate a complete set of eigenfunctions. We would also 
like to emphasize that for the a — > limit, there is a smooth and continuous transition to 
one dimensional harmonic oscillator solutions (even and odd) . The general solutions of eq. 
(|2^ ) in the domain < < 7r/3 are given by 



j2i{(P) ~ COS — sm — F 



2 J V 2 

/2m (0) ~ (sin30)^+'^F 



21,21 + 1,1 - a;cos^ — 



2 

21 - 1,2/ + 2 + 2a,l + a;cos" ' 



2 ^ 
2 



b2i = m + i) , (30) 
3 

b2i+i = 3(2/ + a + -) , 



where a = . The states f2i{4') {f2i+i{4>)) ^'^^ symmetric (anti-symmetric) about the 

point (j) = ^. To determine wave functions for any value of (j), one uses fi{4>) = fi{4> — 
where m is an appropriately chosen integer such that < (p— ^ < ^ . The eigenvalues are 
given by En,2i = \J\{2n + Ql + |)w and En,2i+i = -^/f (2n + 6/ + 3a + ^)uj. 

Discussion and Comments: 

In this paper, we have shown that with a change of parameters of the type a2 = f'^{ao) = 
ao, one gets interesting shape invariant potentials with two shifted sets of equally spaced 
eigenvalues. These potentials are similar to the ones discussed by Calogero|8[]. The main 
differences are: (a) the coefficient a of the inverse square term is always of the "soft" type 
i.e. < a < I and (b) there is an additional 5-function singularity at the origin (in the 
multi-particle case this is a two-body contact interaction). For multi-particle systems that 
interact via two-body potentials of the above type, the spectrum includes the entire set 
eigenvalues of the Calogero potential as well as another set of equally spaced eigenvalues 
which are slightly shifted from Calogero set. 

Can the same technique be extended, for example, to 03 = /^(ao) = ao, thus giving three 
shifted sets of equally spaced eigenvalues? This does not seem to be the case. Not only 
are there no simple real solutions to the constraint /'^(ao) = aQ, but the shape invariance 
condition equations analogous to eq. (^) are not easily solvable. Other choices for the 
change of parameters ai = /(ao) are currently being investigated. 

A.G. acknowledges the hospitality of the UIC Department of Physics where part of this 
work was done. We also would like to thank R. Dutt and A. Khare for many valuable 
discussions and comments. Partial financial support from the U.S. Department of Energy 
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Figure Captions 



Figure 1: The potentials V-{x,ao,e) given by eq. ( |T^ ) for the value e = 0.1 and for two 
choices of parameters: (a) loq = 1.7, toi = 0.3 ; (b) ujq = 0.3, tui = 1.7 . The potentials 
correspond to the superpotentials of Fig. 1. The eigenvalue spectrum, consisting of two 
shifted sets of equally spaced energy levels, is given by eq. (p^) and corresponds to the limit 
e ^ 0. 

Figure 2: Graphical representation of various modes of changing parameters a„ = /"(oo) 
in the shape invariance condition. The variation of a„ versus n is shown by the symbol 
+ for translations a„ = oq + and by the symbol * for scaling a„ = q^UQ. An example 
of the new change of parameters discussed in this paper is /(ao) = 2^ log[coth(Cao)], 
and the results are shown by the symbol o. The numerical choices for this figure are 
ao = 0.8, /3 = 0.1, q = 0.9, C = 0.5 . 

Figure 3: The superpotentials W{x,ao,e) given by eq. (p]). The curves are drawn for a 
small value e = 0.1, in order to see the behavior near the origin. There are two choices of 
parameters: (a) loq = 1.7, ui = 0.3 ; (b) ujq = 0.3, uji = 1.7 . 

Figure 4: The low lying eigenfunctions given by eq. (^) corresponding to the potentials 
of Fig. 1. 

Figure 5: A plot of the five lowest energy eigenvalues Eo to E4 as a function of a = 
{ujQ — u;i)(3u;o + <^i)/4(l<jo + 1^1)^ holding u = {ujq + uii)/2 fixed, a is the coefficient of the 
singular term in eq. (|lO|). 2u> is the equal energy spacing in the two sets of energy levels. 
The two sets are displaced by energy ujq. Regions shown by dashed lines correspond to 
singular but normalizable wave functions. 
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